Mon. Not. R. Astron. Soc. 000,[H15](2002) Printed 1 February 2008 (MN IATpX style file v2.2) 



Camouflaged galactic CMB foregrounds: total and polarized 
contributions of the kinetic Sunyaev Zeldovich effect 



Andre Waelkens 1 , Matteo Maturi 1 ' 2 , Torsten EnBlin 1 

[ Max Planck Institut fur Astrophysik, Garching 85741, Germany 
2 Institut fur Theoretische Astrophysik, Heidelberg 69120, Germany 



o 
o 

(N 



43 

Oh 
6 



o 

O 

o 



x 



Accepted 20?? December ??. Received 20?? December ??; in original form 20?? October ' 



ABSTRACT 

We consider the role of the galactic kinetic Sunyaev Zeldovich (SZ) effect as a CMB fore- 
ground. While the galactic thermal Sunyaev Zeldovich effect has previously been studied and 
discarded as a potential CMB foreground, we find that the kinetic SZ effect is dominant in 
the galactic case. We analyse the detectability of the kinetic SZ effect by means of an opti- 
mally matched filter technique applied to a simulation of an ideal observation. We obtain no 
detection, getting a S/N ratio of 0.1, thereby demonstrating that the kinetic SZ effect can also 
safely be ignored as a CMB foreground. However we provide maps of the expected signal for 
inclusion in future high precision data processing. Furthermore, we rule out the significant 
contamination of the polarised CMB signal by second scattering of galactic kinetic Sunyaev- 
Zeldovich photons, since we show that the scattering of the CMB quadrupole photons by 
galactic electrons is a stronger effect than the Sunyaev Zeldovich second scattering, and has 
already been shown to produce no significant polarised contamination. We confirm the latter 
assessment also by means of an optimally matched filter. 

Key words: scattering - cosmic microwave background - polarization. 



1 INTRODUCTION 

The role of the Sunyaev Zeldovich effect (hereafter SZ) as a CMB 
foreground has been widely studied for c l usters of galaxies (see 
for ex a mple ISunvaev & Zeldovichl i 19721). ISunvaev & Zeldovichl 
1980), iBirkinshawl Jl999h . iDolag et alj jlOOj ) and lSchafer et al.1 
2006)). Due to the high cluster temperatures the relevant effect 



in these cases is the thermal SZ effect (hereafter tSZ), while the 
kinetic SZ effect (hereafter kSZ) is secondary. In this work we 
study the galactic SZ effect. The relevance of the tSZ effect as a 
galactic foreground h as been ruled out by previous studies (see e.g. 
iHinshaw et alj|2007l also section [2.1. It . As we will show, how- 
ever, the kSZ effect makes a significantly higher contribution. This 
is confirmed by an independent g alactic kinetic SZ effect simula- 
tion made by lHaiian et al. I ( l2007h . with whose results ours agree. 
Even though the galactic kSZ effect is more significant than the 
tSZ effect, we demonstrate by means of the optimally matched fil- 
ter technique that it is still not a significant CMB foreground. 

We extend our search by also investigating whether the sec- 
ond scattering of galactic kSZ effect photons could present a sig- 
nificant CMB polarisation foreground. The strength of the the po- 
larised signal due to Thomson scattering by galactic electrons is 
determined by the intensity of the incident quadrupole radiation in 
the rest frame of the electron. Hence we assess the importance of 
a polarised foreground by the intensity of the incident quadrupole 
radiation source. 



We will show, that the polarisation signal due to the intrinsic 
CMB quadrupole is far more significant than that caused by sec- 
ondary scattering of the kSZ induced quadrupole. In addition we 
also examined the relevance of the CMB quadrupole induced by 
the mildly relativistic motion of the galaxy in the CMB rest-frame. 
This effect, however, is known to b e an order of magnitu de smaller 
than the intrinsic CMB quadrupole (Hinsha w et al.l2007l) hence, al- 
though it is slightly more relevant for polarisation production than 
the kSZ quadrupole, also negligible. The CMB quadrupole induced 
polarisation by galactic electrons has been shown to be an insignif- 
icant CMB polarisation foreground dHirata et alj|2005b . We con- 
firm this result in sec. 14.21 where we demonstrate by means of a 
matched filter that no galactic signal can be detected even assum- 
ing an idealised experiment. Hence polarisation contributions due 
to the weaker galactic kSZ effect can safely be ignored as CMB 
polarisation foreground. 



The structure of the paper is the following. In section [2711 we 
give a brief description of the SZ effect, while in section [2~2l we de- 
scribe the Thomson scattering theory and check for the relevance of 
the kSZ effect polarisation signal. Section[3]describes our simula- 
tion of the galactic SZ effect and the polarised Thomson scattering 
emission. In Section[4]we describe our attempt to measure the SZ 
and polarised Thomson scattering signature in our simulations. We 
draw conclusions in Section [5] and discuss how both the kSZ and 
the polarised Thomson signal can be ignored as CMB foregrounds. 
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2 THE SUNYAEV-ZELDOVICH EFFECT 

Here we present a brief theoretical introduction to the SZ effect and 
perform an order-of-magnitude estimate of the galactic tSZ and the 
kSZ effects as a first step toward investigating their importance as 
CMB foregrounds. We do also check whether any relevant contami- 
nation of the CMB polarisation signal due to the galactic kSZ effect 
is expected. 



2.1 Basics 

The SZ effect is a distortion of the CMB signal caused by the 
scattering of the CMB photons on moving free electrons (see 
ISunvaev & Zeldovicd[l972l . 1 198Ch . There are two components of 
the SZ effect: the thermal SZ effect due to the electron thermal 
velocities and the kinetic SZ effect caused by the bulk motion of 
electrons in the CMB rest frame. The CMB spectrum distorted by 
the Compton scattering is given by 



I(x) = Icmb(x) + 8I,sz{x) + SI ksz {x) 
= k U(x) + g(x)y - h(x)rfi] , 



(1) 



where x = hv/(kT) is the dimensionless frequency. The CMB 
blackbody spectrum is 



Icmb(x) = f'o i(x), 



(2) 



with normalisation factor i = 2(kT cmb ) 3 /(he) 2 and spectral shape 
i(x) = x 3 /(e v - 1). The second term on the right hand side of Eq.Q] 
describes the tSZ effect. This effect has a spectral shape given by 



g(x) 



(e x - I) 1 \ e x -I 



e x + 1 
x 4 



(3) 



and amplitude given by the Comptonization parameter for a given 
line-of-sight (LOS), 



y ■ 



m e c l J 



dln e (l)kT e (l). 



(4) 



The last term on the right hand side of Eq. Q] represents the kSZ 
effect. It has the spectral shape 



h(x) ■■ 



(e* - l) 2 
and depends on the LOS integral 



J din, 



mm- 



(5) 



(6) 



Here fj is the velocity component along the LOS in units of the 
speed of light. 

It is common practice to display CMB measurements in terms 
of temperature distortion maps. This is specially convenient for the 
kSZ effect, which can be interpreted as a frequency independent 
change to the original CMB temperature. The spectral distortions 
introduced by the kSZ effect are 

I(x') ^ I(x) + (x' - x) = i - h(x) ^-J , (7) 

and therefore identical to spectral distortions generated by a tem- 
perature change, to first order. Here x' = jp, where T = Tcmb+ST 
and ST <k T C mb- Comparing Eq.Q] and [7] suggest the identification 
Ty6 = ST ksz /T CMB . In contrast, the tSZ spectral distortions cannot be 
described by a simple frequency-independent temperature change. 
However, for our galaxy the kSZ effect dominates by orders of mag- 
nitude over the tSZ effect, as we show in the following. 



2.1.1 Order of magnitude estimate of the galactic thermal SZ 
effect 

Previous works have ruled out the rele vance of the thermal SZ ef- 
fect as a foreground for the CMB (see iHinshaw e t al. 2007). They 
estimate y - 2 x 10~ 8 by assuming a maximal temperature, den- 
sity and LOS distance in Eq. |4j Thus they argue that the SZ effect 
can safely be ignored as a diffuse contaminating foreground signal. 
Indeed the galactic tSZ effect would produce a variation in tempera- 
ture about three orders of magnitude smaller than the CMB primary 
anisotropies, at most. 



2.1.2 Order of magnitude estimate of the galactic kinematic SZ 
effect 

We will now attempt an estimate on the kSZ signal strength. The 
variation in temperature due to the kSZ (in the CMB rest-frame) is 
given by 6T/Tcmb = ?P which is a function of the thermal elec- 
tron density and the LOS components of the electron velocities in 
the CMB rest-frame (Eq. [6^. As a rough estimate, we consider a 
homogeneous gas dist ribution with electron density n e ~ 0.1 cm~ 3 
(Cordes & Lazio 2002), and a constant velocity of 371km/s (see 
lFixsenetal.lll996l) all along the LOS, with a length L = 20 kpc, 
through the galactic plane. Hence, 



T/3 = ST ksz /T CMB ~o- T Ln e j3~4- lO" 6 . 

This toy scenario portrays how large the kSZ effect amplitude 
could maximally be in the direction of the galactic plane. Though 
still smaller than the CMB anisotropy signal, it is significantly 
greater than the thermal SZ effect (about two orders of magnitude). 
This motivates us to undertake a more precise analysis of this phe- 
nomenon by constructing a synthetic map of the galactic kSZ effect 
(see sec. l3.2t and to check whether it can be a non-negligible CMB 
foreground (see sec.|4j- 



2.2 Polarisation due to Thomson scattering 

Polarisation due to Thomson scattering depends only on the 
quadrupole of the incoming radiation as is briefly demonstrated in 
the following. 

The Q and U Stokes parameters due to Thomson scattering 
for each LOS are computed for a local coordinate system whose 
xy-plane is orthogonal to the LOS and whose z-direction points 
to the observer. To consistently compute the Stokes parameters for 
all directions on the sky we adopt the convention where the x-axis 
points in negative 0' and the y-axis in negative & direction of the 
observer-centric coordinate system whose x' axis points towards 
the galactic centre while the z' axis points towards the galactic north 
pole. Assuming the usual 6 and <p angles for the local coordinate 
system we have 



g=^I J dln e (l) J dQ. sin 2 6» cos 20 /(6», 0, /), 
(/ = ^I J dln e (l) J dQ. sin 2 6> sin 20 /(#, 0, /). 



(8) 
(9) 



The terms sin 2 9 cos 20 = V8*715 (y, 2 (6», 0) + Y~ 2 (9, 0)) and 

sin 2 6» sin 20 = -i V87r/15 (y 2 (9, 0) - Y; 2 (6,(pj} are linear combi- 
nations of 1=2 spherical harmonics, i.e. the quadrupole (here spher- 
ical harmonics are denoted by Yj"(9, 0)). Hence only the quadrupole 
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part of 1(8, (p) survives the integration over the whole sphere: 



f 



dQ. sin" 8 cos 20/(9,0) 



and similarly 

J dQ sin 2 8 sin 741(6, 




(a 22 -a 2 _ 2 ) 




Im a 22 . 



(10) 



(11) 



Here, the a\ m are the spherical harmonic coefficients. We have made 
use of the property fl;_,„ = —V"a* lm of the real valued brightness 
distribution, with the asterix denoting complex conjugation. 



Quadrupole source 



C 2 



CMB 

relat. Doppler effect 
kSZ effect 



221^K 2 
3.45^K 2 
l(T 2 /iK 2 



Table 1. Amplitude of the quadrupole moments in order of their signifi- 
cance. The relativistic Doppler effect quadrupole moment is cal culated for 
j8 ~ 1.2 ■ l(r 3 . The CMB quadrupole given by iHinshaw et all fe007l) in- 
cludes this value (i.e. they do not attempt to subtract it from the measured 
value). The kSZ quadrupole moment is obtained from our simulation in sec. 
13.21 Other foregrounds with frequency dependent brightness temperatures 
like synchrotron or dust emission are assumed to be removed. 



consequently neither is the one due to the second scattering of kSZ 
photons. 



2.2.1 Estimating the relevance of the polarised signal 

Here, we asses the typical strength of the polarisation signal caused 
by second scattering of the kSZ photons. To do so we compare it 
to other possible sources of polarised radiation. In order to identify 
the dominant polarisation signal we only have to find the strongest 
quadrupole source. 

TableQ]lists the total power of all possibly relevant quadrupole 
sources. Notice that other foregrounds like synchrotron or dust 
emission present even stronger quadrupoles, however, since they 
have frequency dependent brightness temperatures, they can, in 
principle, be removed. We see that the CMB quadrupole is by 
far the most significant contribution. The Doppler quadrupole due 
to Earths motion in the CMB rest-frame is negligible, while the 
quadrupole due to second scattering of the kSZ effect is even 
smaller. 

We first perform an order of magnitude estimate to obtain the 
upper limit for the polarised signal due to the scattering of the CMB 
quadrupole by galactic thermal electrons: 

W e refer to the CMB power spectrum given by Hinshaw et al. 
J2007I) . In their table 7, they display the quantity AJf = 21 \pK 2 . It 
is given by AT 2 = /(/ + l)Q/(2n), and C, = 1/(2/ + 1) Z' m= -i l"U 2 - 
For a rough upper limit estimate, we assume all quadrupole power 
to be concentrated on only the real part of a 22 = V5C7. Assuming 
a homogeneous thermal electron num ber density n e ~ 0. 1 ctrr 3 
(again following ICordes & Lazidr2002l as in section r2. 1.2t along 
a LOS of L ~ 20 kpc (say towards the galactic centre) we get, 



/ Ln e Kea 22 m 2 ■ 10 _ pK. 



(12) 



The polarised CMB intrinsic anisotropies have maximum ampli- 
tudes of ~ 10% of the total CMB anisotropies whose maxima ar e 
typically around a few 100pK (see figure 15 o f IPage et alj|2007r) . 
Hence, the polarisation signal due to Thomson scattering on the 
galactic electrons from the CMB radiation (only a few 10~ 2 //K) is 
likely to be undetectable on any Q or U maps. The much weaker 
kSZ contribution (see table [T} is thus likely to be an insignif- 
icant contamination to the CMB polarisation. This is supported 
by detailed computations of the Thomson sc attering of the CMB 
quadrupole by galactic electrons performed bv lHirata et all (2005). 
The authors verify that the galactic Thomson signal is not a signif- 
icant contamination to the intrinsic CMB polarisation. 

Complementary to the arguments presented in iHirata et"al] 

d2005h . in sec. |4.2| we confirm by means of the matched filter tech- 
nique applied to our simulations that the signal is not detectable, 



3 SIMULATIONS 
3.1 The code 

We modified the Hammurabi-code origin ally designed for sim- 
ulatin g galactic synchrotron emission jWaelkens et"ai] 120071: 
EnBlin et al. 2006), to calculate the LOS integral given in Eq. [6] 
(for the kSZ effect) and in Eq. [8] and [9] (for the polarisation due 
to Tho mson scattering). Th is code uses the HEALPix pixelisation 
scheme dGorski et ai]|2005l) . 



3.1.1 Line-of-sight integration for the kSZ effect 

The flux associated with an observation beam, whose width and 
shape is related to the HEALPix pixels, is calculated by taking sam- 
ples of the thermal electron density and the velocity along the LOS 
at equal radial intervals. To compensate for the widening of the 
beam cross section as one moves farther away from the observer, 
the observation beam is, after a certain distance, split into 4 sub- 
beams over which we then average. Each sub-beam is again indi- 
vidually sampled at the same regular radial intervals as before, and 
split again after twice the first distance, and so forth. This avoids an 
under-sampling of the thermal electron density and LOS velocity 
at large distances from the observer. 

3.1.2 Line-of-sight integration for the polarised galactic 
Thomson scattering 

In this case only the electron density is integrated along the line 
of sight, however the computation is performed in an analogous 
manner to the previous paragraph. The resulting map is multiplied 
by the a 22 component after a rotational coordinate transformation 
for each different LOS in order to get the respective Q (Eq.[8} and 
U (Eq. [9} Stokes parameters. The computation of the proper a 22 
coefficients is described in appendix 151 



3.1.3 Galaxy model 

Two input models are needed, a model for the thermal electron dis- 
tribution and a model for the galactic rotation curve. 

(i) For the elect rons we use the NE2001 model (see 
ICordes & Laziol2002l) . This model subdivides the galaxy in several 
large scale structure elements (thin disk, thick disk, spiral arms,...) 
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Figure 1. ST = rfi ■ T map in /jK. The monopole and dipole components of 
the map have been subtracted in analogy to what typically is done for the 
observations (e.g. WMAP). Our simulation shows that the maximum ampli- 
tude of the map decreases by about 20% once the kSZ dipole is subtracted. 
Notice the features due to over dense or under dense regions in the local 
interstellar medium. 



Figure 2. Stokes Q in /iK. 




as well as some local small scale elements such as supemovae bub- 
bles. 

(ii) For the galactic rotation curve we assume a constant value 
of 220km/s above a galactocentric distance of 2kpc. Below this 
distance the velocity falls off linearly to zero while approaching the 
galactic centre. This is an approx imation for the observed rotation 
curve (see e.g. iKlvpin et alJl200l and references therein), which 
should be adequate to determine whether the galactic kSZ effect 
can be detected. 



3.2 Constructing the kSZ map 

The velocity of each point of the galaxy with respect to the CMB 
rest-frame is subdivided into two components. The rotation veloc- 
ity of the galaxy (p, ot ) and the position-independent translation ve- 
locity of the galaxy (j8„„). The first quantity is obtained from our 
rotation curve model. The later quantity is obtained by subtract- 
ing our solar systems galactic rotation velocity vector (220 km/ s in 
/ = nil direction) from the measured total CMB fr ame velocity 
(371k m/s in / = 1.477T and b = 0.27tt direction, see iFixsen et alj 
1996). Since the resulting /? is small(~ 1CT 3 ), we can safely ignore 
relativistic effects when adding velocities. 

However, relativistic effects may need to be considered for the 
temperature transformation. So far we obtained the t/3 = 6T/Tcmb 
sky in the CMB rest-frame. The relativistic temperature transfor- 
mation is given by 8T obsf = 6T cmbf \y{\ + (j8 to ,jU oi , s/ ))] . Here 
Hobsf i s tne cosine of the angle between the opposite of the CMB 
frame motion direction and the LOS in the observer rest-frame and 
7 = (1 - /3 2 y >12 . The relativistic correction to the temperature is 
small, in the extreme a factor ~ (1 ± 10~ 3 ), which can be safely 
ignored, leaving 6T obs f = ST C mb} - Thus to obtain the kSZ map we 
can neglect all relativistic effects. 

Finally, we obtain the kSZ map by inserting our galaxy model 
in Eq. [6] The resulting map is shown in figure [T] Note the local 
interstellar medium features of the thermal electron density model 
(see sec. 13.1.31 . We subtracted the monopole and dipole in accor- 
dance with the standard practice of subtracting the monopole and 
the relativistic-Doppler-effect induced dipole from maps displaying 
the CMB anisotropies. Our simulations show that this has a small 
effect, decreasing the maximum amplitude of the map by about 
20%. 

The maximum amplitudes of the simulated signal are of the 




Figure 3. Stokes U in fiK. 



order of pK, in accordance with the estimate obtained in section 
12.1.21 and two orders of magnitude smaller than the CMB intrinsic 
anisotropies. 



3.3 Constructing the polarised Thomson scattering maps 

The expressions for the Q and U stokes parameters (Eqs. [8] and 
[9]l assume that the negative z-axis of the coordinate system co- 
incides with the LOS direction. Thus it is necessary to compute 
the appropriate a 2 2 coefficient for each direction on the sky. This 
is done by rotating a set of o^m coefficients, which are given in 
a certain coordinate system e.g. that in which the WMAP data 
are presented, to the proper coordinate system in which the neg- 
ative z-axis is aligned with the LOS. The procedure for rotatin; 
a hn coefficients is described in det ail in lVarshalovich et"aT] dl98 
and de Oliveira-Costa et al. I d2004h . For completeness, we repro- 
duce our special case in the appendix [B] Having performed the ro- 
tation we may apply Eqs. 1 101 and II II to each position on the sky, 
which after performing a LOS integral of the thermal electron den- 
sity, finally allow us to compute Q and U, using Eqs. [8] and [9] re- 
spectively. The resulting full sky maps are shown in figures [2] and 
[3] The maximum amplitudes are ~ 10~ 3 yi/K, about one order of 
magnitude lower than our estimate in section l2~2l Out of the Q and 
U maps, the E and B mode spherical harmonic coefficients are ob- 
tained. Since most of the CMB polarisation power is in E modes, 
the B-modes signal of the polarised galactic Thomson scattering is 
more relevant as a CMB polarisation contamination than the corre- 
sponding E-mode signal. In sec. 14.21 we search for B-mode signal 
contamination by polarised galactic Thomson-scattering. 
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Figure 5. Simulated CMB sky contaminated by the simulated galactic kSZ 
effect. 



Figure 4. Comparing the CMB power spectrum to the simulated kSZ power 
spectrum. The saw like pattern for the kSZ effect power spectrum is due to 
the fact that even spherical harmonics are zero at the equator of the coor- 
dinate system, which in this case corresponds to the galactic plane. Hence 
odd modes will have more power. 




4 FILTERING 

In order to estimate the feasibility of a detection of the effects dis- 
cussed above, we construct and apply an optimally matched filter 
to simulated data, where we probe for the galactic kSZ and galactic 
polarised Thomson scattering signatures. 

The optimal matched filter was proposed for measuring ra- 
dial bulk motion of galax y clusters through their kSZ imprint 
faaehnelt & Tegmarki 1996) and subsequently found a wide range 
of applicatio ns in the literature, such as weak lensing and the RS ef- 
fect (see e.g. iMaturi et all2 005. 2006). A sc alar adaptive generali - 
sation of the filter on the sphere was done bv lSchafer et al j (2006), 
and later ext ended to allow non-a zimuthally symmetric objects to 
be detected ( McEwen et al ] l2006h . We use a similar approach to 



ISchafer et al 



(2006). However, instead of convolving the CMB sig- 



nal and the filter for various positions and orientations, we need 
only perform one convolution since we already know the position 
and orientation of our signal. Details are presented in the appendix 
lAl The estimate of the maximum amplitude of the targeted signal 
is given by the convolution of the filter T with the signal 5 



while the variance is 
a- 2 = <A> 2 



(13) 



(14) 



Here \j/i m is the spherical harmonic coefficient of the filter. The an- 
gular power spectrum of the noise (C/) is either, in the kSZ case, the 
temperature angular power of the CMB plus instrumental Gaussian 
noise, or, in the polarised Thomson-scattering case, the B-mode 
angular power spectrum plus corresponding Gaussian instrumental 
noise. The filter is defined such that a 2 is minimal while the bias 
b = (A) - A = is zero. 

4.1 Filtering on simulated data: galactic kSZ 

We consider the ideal situation in which the only additional contri- 
bution to the CMB plus instrumental noise signal is the kSZ effect. 

We do a simulation of the CMB sky with a sampling resolu- 
tion of about 7 arcmin (corresponding to the HEALPix parameter 




Figure 6. The optimally matched galactic kSZ filter. Note the features due 
to the l ocal inter stellar me dium components of the thermal electron density 
model ICordes & Laziol2002h . 



NSIDE=512) and assume experimental characteristics expected for 
a Planck-like experiment, i.e. a resolution of 10 arcmin and a sen- 
sitivity of 5yuK/beam. Given that the structures of the simulated 
kSZ effect are mostly on larger scales (see figure[4]l, the sensitivity, 
instead of the resolution, is the more relevant instrumental charac- 
teristic here. We generate a kSZ map with the same sampling reso- 
lution, add it to the synthetic CMB sky and then convolve the map 
with a one degree Gaussian kernel. This smooths the sharp borders 
of the simulation and also ensures that our final map is properly 
sampled. The result is shown in figure[5] in which no evidence of 
the kSZ effect is visible. 

As a signal template we take the same convolved kSZ map 
previously used in defining the data. This highly idealised approach 
(assuming a perfect signal template and no contamination by other 
Galactic emission sources) is suitable to determine the best possible 
filter performance. 

The signal amplitude obtained by applying the filtering tech- 
nique is -10 ± 21//K, with the variance given by Eq.[l4] The ex- 
pected maximum amplitude of the kSZ map, the quantity we are 
estimating with the optimally matched filter, is 2.4/iK, according to 
our model. Hence we have a signal-to-noise ratio of 0.1. 

Although in this simulation no other foregrounds were 
present, no detection was possible, demonstrating that the galactic 
kSZ effect can be ignored as a CMB foreground. 



© 2002 RAS, MNRAS 000,[[H9] 



6 Andre Waelkens, Matteo Maturi & Torsten Enfilin 



4.2 Filtering on simulated data: the B-mode signal 

The measurement of CMB B modes is the next cosmological chal- 
lenge and would be very relevant for constraining inflationary mod- 
els. The CMB polarised power is overwhelmingly in the E-modes, 
contrary the polarised galactic Thomson scattering. For this, our 
simulation shows a balance between the E and B-mode power on 
large scales, and a difference by at most one order of magnitude on 
smaller scales. This implies that the galactic B-mode contamination 
of the CMB is relatively stronger than the corresponding E-mode 
contamination. Therefore applying the same filtering technique as 
in sec. [4] we can investigate whether the polarised Thomson scat- 
tered B-mode signal is detectable and, if so, whether it is a relevant 
CMB foreground. 

We use the predicted CMB BB power spectrum Cf B (as given 
by the CAMB code. 



Lewis et al. 



2000) for a tensor-to-scalar ratio 



of 10 5 (following iHirata et alj|2005h . plus the instrumental noise 
for a sensitivity of 7.6 /jK/beam and a resolution of 30 arcmin (cor- 
responding to polarisation instrumental characteristics of a Planck- 
like experiment), as the noise power spectrum in Eq. |A10| The ten- 
sor spherical harmonic coefficients — > j Un for the B modes 
of our polarised galactic Thomson-scattering simulation are used 
as the spherical harmonic coefficients of our template. 

This analogy to the previous case is permissible since Cf B 
is completely uncorrelated between each different mode (like the 
CMB temperature power spectrum), and the CMB B mode signal 
is uncorrelated with the galactic B mode signal. 

We obtain an estimate of -4 • 10~ 3 ± 5.8 • 10~ 2 yuK, against 
the simulated signal amplitude of 1.6 • 10~ 3 /iK. Since the achieved 
signal-to-noise ratio is ~ 0.03, the Thomson scattering by galactic 
thermal electrons is a negligible contribution to the CMB polarisa- 
tion signal. 



a tensor to scalar ratio of 10~ 5 (where most of the polarised CMB 
power is in E-modes, and hence the B-modes are more easily con- 
tamined). The polarised galactic Thomson-scattering signal is also 
concentrated on the galactic plane. The maximum amplitude of the 
simulated signal is 8 • 10~ 4 /iK. Assuming a Planck-like instrument 
with a resolution of 30 arcmin and a sensitivity of 7.6 yuK/beam, we 
estimate an amplitude of -4 • 10~ 3 ± 5.8 • 10 _2 yuK; thus no contam- 
ination of the CMB B-mode signal is detectable. 

In conclusion, galactic foregrounds due to the kinetic 
Sunyaev-Zeldovich effect are a certainly existing, however per- 
fectly camouflaged contamination of the CMB. Nevertheless, it 
is an expected physical effect which generates 1% temperature- 
fluctuation amplitude corrections on large scales, and thus should 
be taken into account when estimating the error of precision mea- 
surements as the Planck satellite will provide. 
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5 CONCLUSIONS 

We have investigated the role of the galactic kinetic Sunyaev- 
Zeldovich effect as a CMB foreground, both for the total and po- 
larised emission. This was done by means of an optimally matched 
filter based on a model of the signal and the expected noise power 
spectrum. The filter returns an unbiased estimate of the peak signal 
amplitude and is optimal in the sense that it provides the minimum 
variance for such a measure. 

Regarding the non-polarised kinetic Sunyaev-Zeldovich emis- 
sion as a CMB foreground, we found that, although the ki- 
netic Sunyaev-Zeldovich effect is more relevant than the thermal 
Sunyaev-Zeldovich effect, it is still not strong enough to be de- 
tectable by means of an optimally matched filter. 

Our simulation for the kinetic Sunyaev-Zeldovich effect 
shows the signal to be concentrated on the galactic plane where 
the electron density is highest. The maximum amplitude of the 
simulated signal (the quantity estimated by the filtering procedure) 
is 2.4/jK. Assuming an instrument with resolution of 10 arcmin 
and a sensitivity of 5 /jK/beam, we obtain a signal amplitude of 
-10± 21/jK. Hence, no kinetic Sunyaev-Zeldovich effect signature 
could be detected. 

We further demonstrate that the secondary scattering of kinetic 
Sunyaev-Zeldovich photons is negligible as a polarisation source. It 
is a weaker effect than the polarised emission due to Thomson scat- 
tering of CMB photons by galactic electrons. Again by means of the 
optimally matched filter technique, we determined that the latter ef- 
fect is also irrelevant even as a CMB B-mode foreground, assuming 
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APPENDIX A: FILTER DERIVATION 

In this section we describe how to obtain an optimally matched 
filter function ¥(0, (p) = Z m ^imY'"(6, (j>) for attempting to detect 
a non azimuthally-symmetric target signal (such as the signal due 
to e.g. the galactic kSZ or the B mode due to galactic Thomson 
scattering) in the presence of a Gaussian random noise (e.g. the 
sum of CMB anisotropics, total or polarised, and the instrumental 
noise). W e used such a filter ing procedure in investigations of the 
RS effect dMaturi et alj|2007l) . 

F o llowing the argume nts l ayed out by Haehnelt & Tegmarkl 
dl996l) , iMaturi et alj |2003) and ISchiifer et al.l J2006I) . we assume 
the data to be composed of the statistically homogeneous CMB 
anisotropics, some instrumental Gaussian random noise and the 
targeted signal. The expansion of the data in spherical harmonics 
(17(0,0) ^ 



; m 

where we model the complex coefficient j/„, as 
$lm = At/,,, + a lm . 



(Al) 



(A2) 



We choose A to be the maximum amplitude of the targeted sig- 
nal, while the morphology is encoded in the complex spherical har- 
monic coefficient r lm . The sum of the CMB and the noise sky is 
ay, = a 'i,„ + n im, where «/,„ is the instrumental noise modelled as a 
white Gaussian random noise field. Note that the ensemble average 
(aim) = 0. 

In the studies mentioned above, the convolution of the filter 
with the total signal is performed for all positions on the sky (or 
the corresponding analogous on the plane), in effect generating a 
convolution map described by 



¥*S= da*¥(0-a,<l>-P)S{0,(l>). 



(A3) 



In our case, however, we know the precise position and orientation 
of our target signal, therefore we need only to evaluate the convo- 
lution for a single position on the sky. We define the estimate of the 
target signal strength to be 



dns(6, <pyp(o, (/>). 



(A4) 



We want to estimate the real peak amplitude of the target effect 
accurately. Hence we enforce an unbiased estimate 



b = A- < A >= 0, 
with variance 

-2 



a 1 = (A )-A 2 = Y,Yj C '^ 2 



(A5) 



(A6) 



We enforce the condition that the variance be minimal by means of 
Lagrange multipliers. We are looking for the filter Y which would 
give the smallest error for the estimate in Eq. IA4I subject to the 
constraint of being unbiased (Eq. lA5t . Thus we solve 

-JL(<r 2 -Ab-A'b-)=0. 

The resulting filter is given by 

Ati„, 



(A7) 



tylm 



Ci 



which, when substituted in Eo. lA5l yields 
1 



A = 



ZiZ,„{\t,J 2 /c) 



and 



(A8) 



(A9) 



complex value 



a 2 2 (-14.41 -il8.8)/V2 
a 2l (-0.05 + (4.86)/ VI 
a 20 11.48 



Table Bl. T he a/ m coefficie nts calculated from the real valued coefficients 
in table 6 of Hins haw et al l <2007t) . These coefficients are obtained for the 
z-axis of the coordinate system pointing to the galactic north pole. 



(A10) 



Z,;E„,( |T "» |2/C ') Cl ' 

Analogous to the flat space results in lHaehnelt & Tegmarkl Jl99rj) . 
the filter gives preference to modes where the ratio of the target 
signal to the noise is largest. The estimate of the target signal am- 
plitude is performed using Eq. lA4l with the error given by IA6l and 
the filter bv lATOl 



APPENDIX B: ROTATING SPHERICAL HARMONIC 
COEFFICIENTS 

Here we describe how a set of spherical harmonic quadrupole coef- 
ficients (a 2m ) transforms under a rotation of the coordinate system. 
In particular we obtain the rotated spherical harmonic coefficient 
needed for computing the Stokes Q and U parameters (Eqs.[8]and 
|9]in sec.im. 

The real-valued a 2m coefficients for the CMB radiation given 
bv lHinshaw et alj d2007l) , are related to the complex-valued spher- 
ical harmonics (the convention in this paper) as listed in table |B1 

The transformation of spherical harmonics under a rota- 
tion described by the E uler angles a, /3 and y is given by 
IVarshalovich et al.ldl98^) : 



(Bl) 



where D' mm ,(a,[}, y) is the Wigner-D symbol. The hat denotes quan- 
tities in the rotated coordinate system, otherwise quantities are 
taken to be in the usual coordinate system where z points towards 
the galactic north pole. In our case, the inverse transformation is 
needed, 



Yf(9,4,) = ^»;„^(a,# r )y;"(^). 

m=-Z 

The inverse of the Wigner-D symbol is given by 

D- ll mm ,(a,j3,y) = D m , m \ a ,/3,y), 



(B2) 



(B3) 



here the asterix denotes complex conjugation. 

We can write the total quadrupole intensity as a combination 
of spherical harmonics in any coordinate system 



//= 2 (0, 0) = J] aim Y ™ (6 ' ® = Z him 



(B4) 



Substituting IB2I in IB4I multiplying by Y 2 (9,(f>), and perform- 
ing an integral over the whole sphere yields, 



"22= ^ a 2m >Dl, 2 "(a,f3,y). 

m'=—2 



(B5) 
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The Wigner-D symbol is D 2 m , 2 *(a,f3, y) = e im ' a dl, 2 *(fi)e i my, with 

* Vm rno. 'M« 'MAH'/a cos <2W >Q8/2) S in (2 -'>(/3/2) 

4„'2 03)= (2 + m )!(2-m )!4! . 

(2 - m )!(2 + m')! 

We require the value of S22 in a coordinate system which has 
its negative z-axis pointing in the direction specified by the usual 6 
and angles. The corresponding Euler angles are 

7 = -7r/2 , J3 = n — 6 , a = (/>- 7T. 

Note that this means that the x-axis and y-axis point in negative (f> 
and negative directions, respectively. This is the coordinate sys- 
tem in which the Stokes Q and U parameters are calculated in this 
work. 

Inserting I122 in Eqs.[lO]and[TT|and multiplying these with the 
corresponding LOS integral of the thermal-electron density and the 
appropriate constant pre-factor yields Eqs.[8]and[9] 

This paper has been typeset from a Tj5X/ KTjSf file prepared by the 
author. 
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